Driven optical lattices as strong-field simulators 



Stephan Arlinghaus and Martin Holthaus 
Institut fiir Physik, Carl von Ossietzky Universitdt, D-26111 Oldenburg, Germany 

(Dated: June 21, 2010) 

We argue that ultracold atoms in strongly shaken optical lattices can be subjected to conditions 
similar to those experienced by electrons in laser-irradiated crystalline solids, but without introduc- 
ing secondary polarization effects. As a consequence one can induce nonperturbative multiphoton- 
like resonances due to the mutual penetration of ac-Stark-shifted Bloch bands. These phenomena 
can be detected with a combination of currently available laboratory techniques. 
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I. INTRODUCTION 



The investigation of ultracold atoms in optical lattices 
constitutes a major area of topical research [ll-I3|- One 
of the long-term visions driving this trend stems from 
the prospect of using these well-controllable and flexible 
systems for "emulating" important quantum many-body 
problems which still are not fully understood, such as 
high-Tc superconductivity [5;, ,6], and of obtaining infor- 
mation on these by observing their cold-atom-emulated 
versions in the laboratory, rather than attempting nec- 
essarily imperfect computer simulations. So far, inter- 
est has been focused mainly on systems governed by a 
time-independent Hamiltonian operator, a hallmark ex- 
ample being provided by the Bose- Hubbard model 0- 
However, it is feasible to subject the lattice atoms to 
time-dependent external forces, and thus to study explic- 
itly time-dependent phenomena [1, [!]■ Already in 1998 
Madison et al. have obtained evidence for Bloch band 
narrowing with cold sodium atoms in time-periodically 
forced optical lattices [l^l ; more recently, dynamic local- 
ization [ll|, [l^ , photon-assisted tunneling [l3| , and co- 
herent control of the superfluid-to-Mott insulator tran- 
sition [l3l have been demonstrated with Bose-Einstein 
condensates in such strongly shaken periodic potentials. 
Moreover, it has been suggested to employ oscillating op- 
tical lattices for realizing frustrated quantum antiferro- 
magnetism [l5| . In this article we argue that ultracold 
atoms in forced optical lattices also lend themselves to 
the study of multiphoton-like transitions under strong- 
field conditions which are barely accessible with elec- 
trons in solids irradiated by high-power lasers; in partic- 
ular, they provide an exceptionally clean testing ground 
for the investigation of nonperturbative multiphoton-like 
resonances. We first sketch in Scc.|ll]the required setup, 
and specify the orders of magnitude of the relevant pa- 
rameters which characterize the optical-lattice analogs 
of strong laser fields. We then present numerical model 
calculations in Sec. IHIl demonstrating how both pertur- 
bative and nonperturbative resonances manifest them- 
selves. The explanation of these phenomena makes use of 
both the spatial periodicity of the optical lattice and the 
temporal periodicity of the driving force: Effectively, one 
encounters a spatiotemporal crystal, the band structure 



of which is controlled by the parameters of the driving 
force. This viewpoint is emphasized in the concluding 
SecHVl 



II. SIMULATING STRONG LASER FIELDS 

A one-dimensional (ID) optical lattice is created, for 
example, by shining laser radiation with wavelength A = 
27r/fcL against a mirror and retroreflecting the beam into 
itself. An atom of mass M moving in this standing 
light wave then experiences a periodic potential with a 
depth Vb which is proportional to the laser intensity [1]. 
Mounting the mirror on a piezoelectric actuator now al- 
lows one to let it oscillate sinusoidally with a precisely 
controlled angular frequency uj and amplitude L, thus 
shaking the lattice back and forth fl3| . In the laboratory 
frame, the Hamiltonian describing the particle's center- 
of-mass motion along the lattice direction then reads 



-f^iab = ^ + Y cos {2fcL [x 



Lcos{ujt)]} . (1) 



The relevant characteristic energy scale is given by the 
single-photon recoil energy. 
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typical scaled lattice depths Vo/E^ range between about 
5 and 10. For example, with ®''Rb atoms in a lattice 
erected by light with wavelength A = 842 nm one has 
E^ — 1.34 X 10"^"'^ eV, as corresponding to the recoil 
frequency — E,-/{2Trh) = 3.23 kHz. 

Performing a unitary transformation to a frame co- 
moving with the lattice, the Hamiltonian acquires the 
suggestive form [l3. El] 
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cos{2ki^x) — Fx cos{ujt) 



(3) 



with F = MLuP denoting the amplitude of the inertial 
force appearing in this oscillating frame. A meaningful 
measure for the strength of this force is the dimensionless 
parameter [l2j 
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FIG. 1: Lowest three Bloch bands of a ID optical lattice with 
depth Vo — 5.7 Er. The lowest band gap is 2.763 iJr at fc = /cl. 

where d = X/2 specifies the lattice constant. In terms of 
quantities directly accessible in the laboratory, one has 

TT^ v L 

Ko = —--, 5 
2 i^i- a 

with the driving frequency ly — uj/{2tt), showing that 
one may easily realize values Kq > 1 when both ratios 
jz/j/r and L/d are on the order of unity [lll - [l^ . To ap- 
preciate what this means, consider an atomic analog: A 
common KrF exciplex laser provides photons with en- 
ergy hu! — 5.0 eV. Inserting this into the expression Q, 
taking the Bohr radius for the length d, and solving for 
the electric field strength £ — F/e acting on an elec- 
tronic charge, one finds that Kq = 1 is reached only 
for £ = 9.45 X lO^*^ V/m, which is roughly one-fifth of 
the field formally experienced by a ground-state electron 
in the hydrogen atom. In this sense, time-periodically 
forced optical lattices can serve even as superstrong-fleld 
simulators: Shaking a lattice with large amplitudes L 
according to the Hamiltonian ([T]) simulates perfectly ho- 
mogeneous fields acting on particles in periodic poten- 
tials in the regime Kq > 1 of the parameter (|4]) which 
is hard to reach with laser-driven electrons in traditional 
solids, without introducing, for example, detrimental po- 
larization effects. Thus, ultracold atoms in driven optical 
lattices offer the unique possibility to study superstrong- 
field-induced multiphoton-like processes in periodic po- 
tentials in their purest form. 

III. PERTURBATIVE AND 
NONPERTURBATIVE MULTIPHOTON 
TRANSITIONS 

For illustrating the dynamics that become explorable 
in this way, we consider a ID lattice with depth 



Vo = 5.7 El-. Its single-particle eigenstates are Bloch 
waves [17i ]. 

iPn,k{x) = exp(ifca:)u„,fe(x) , (6) 

with lattice-periodic functions Un.k{x) = Un,k{x d) la- 
beled by a band index n and a wave number fc; Fig. [1] 
depicts the energy dispersion relations En{k) for the low- 
est bands n = 1, 2, 3. In the center of the Brillouin zone, 
that is, at k/ki^ = 0, one has E2{Q) - Ei{Q) = 4.690 -Br, 
and i?3(0) - Ei{Q) = 5.544 i?r. We now take an initial 
state exclusively populating the lowest band, as described 

by 

^p{x,tn) = / dkgi{k,to)ipi,k{x,tn) (7) 

J-ki, 

with a Gaussian fc-space distribution 



gi{k,to) = (2fcLV^AA;)-i/2 




centered around k/ki^ — 0, and set Ak = 0.1 /cl for its 
width, as appropriate for an initial ensemble of noninter- 
acting ultracold atoms. This state then is subjected to 
pulsed forcing with an amplitude F(t) which rises from 
zero to a maximum value, stays constant for a while, and 
decreases back to zero. For the sake of definiteness, we 
consider conditions as already realized experimentally in 
Ref. [3] : We take *^Rb as atomic species in a lattice with 
A = 842 nm, and design the envelope of the pulse such 
that F{t) rises linearly within 10 ms, stays constant for 
a holding time of 2 ms, and then is linearly switched off 
in another 10 ms. For a driving frequency of 5 kHz, say, 
the ramp time of 10 ms corresponds to 50 cycles, so that 
the relatively slowly changing envelope F(t) may enable 
adiabatic following under nonresonant conditions. 

Moreover, we rely on the fact that the fraction of 
atoms surviving in the lowest band can be accurately 
determined, as demonstrated by the Landau-Zener mea- 
surements reported in Ref. [18]. We therefore compute 
the escape probability from the lowest band after each 
pulse, for specified values of K^^^ reached during the 
plateau phase. Figure [2] shows results thus obtained for 
i^max ^ Q 7 g^^^ ^max ^ 3^ fuuctious of the driv- 
ing frequency a;/(27r). The pronounced peak pattern de- 
pends markedly on the maximum driving amplitude; for 
instance, a further peak has appeared for K^^^ — 1.3 at 
uo/{2tt) w 4 kHz which was not visible for X™^'' = 0.7. A 
more complete picture is provided by Fig. [31 which shows 
a two-dimensional plot of the escape probability consid- 
ered as function of both w/ (27r) and K^^^, for the same 
pulse shape as taken in Fig. [51 

The positions of most of the peaks in Figs. [H and [3] 
(i.e., most of the system's resonant frequencies) are easily 
explained: Because the initial state is narrowly centered 
around k/ki^ = 0, its response is mainly determined by 
the energies En{0) in the Brillouin-zone center. Hence, 
one expects ordinary m-photon-like resonances between 
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FIG. 2: (Color online) Escape probabilities from the lowest 
Bloch band after pulses with linear switch-on and switch- 
off ramps of 10-ms duration each, and a holding time of 
2 ms, during which a specified value Kq^^^ of the scaled 
amplitude Q is reached. Driving frequencies uu/{2tt) cor- 
respond to *^Rb in an optical lattice with A = 842 nm. 
Light, K^'''' = 0.7; black, Jf^T'"' = 1-3. Of particular in- 
terest is the unexpected, strong, and narrow resonance at 
cj/(27r) = 5.3 kHz. 




FIG. 3: (Color online) Escape probability versus both driving 
frequency ui/{2-k) and maximum scaled amplitude A'™'"', for 
the same pulse shape as employed in Fig. (2] 



the initial band n = 1 and higher bands n — 2,3,... 
when the driving frequency complies with the condition 



A£:„,i = £:„(0) - £:i(0) = mhw 



(9) 



for integer m. Indeed, listing these expected m-photon 
transition frequencies in Table |T] and comparing them to 
the frequencies of the peaks observed in Fig. [31 one gen- 
erally finds quite good agreement. 

In some instances, however, the numerical solution of 
the Schrodinger equation produces a peak which does 
not fit into this naive pattern. Most notably, the sharp 
spike visible in Fig. [5] at uj/{2-k) = 5.3 kHz does not 
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TABLE 1: Expected and computed resonance frequencies: u^cs 
is an m-photon transition frequency according to Eq. ([9}, 
fpcak is the position of the corresponding peak where it be- 
comes apparent in Fig. O The entry NV indicates that no 
peak is visible for the pulse profile employed here. 
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FIG. 4: Quasienergies e„(fc) for the ID optical lattice driven 
with frequency uj/{2tt) — 5.30 kHz, and scaled amplitudes 
Ko = 0.7 (a), 1.0 (b), and 1.3 (c). The areas shaded in 
gray, extending from k = — O.I/cl to A; = +0.1 /cl, mark the 
range of wave numbers explored by the initial wave packet. 
The insets show how the quasienergy band n — 1 (above) 
is pinched through with increasing Kq by the band n = 2, 
displaced downward by Sfku. This causes the nonperturbative 
resonance observed in Fig. (2] 



match Eq. ([9]) for any reasonable combination of n and 
TO. Such "nonperturbative" events are our main concern; 
we predict that they can be detected experimentally in 
already existing setups. These particular resonances ad- 
mit a systematic explanation which forces us to go way 
beyond the perturbative reasoning underlying Eq. 

Because the Hamiltonian ^ is periodic both in space 
(with lattice period d = tt/Ajl = A/2) and in time (with 
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driving period T = 27r/ixi), it gives rise to spatiotemporal 
Bloch waves [11] , 

ipn.k{x, t) = Un^k{x, t) exp {\[kx - en{k)t/h]} , (10) 

with functions Un,k{x, t) = w„,fc(a; + d, i) = Un,k{x, t + T) 
reflecting translational invariance in space and time on 
equal footing, and quasienergies Snik), in generalization 
of the usual Bloch waves ^ . While quasimomenta hk are 
determined up to an integer multiple of 2Trh/d = 2hki^, 
quasienergies are likewise determined up to an integer 
multiple of the photon energy 2-Kh/T = fvjj. Figure |3] 
shows one "quasienergy Brillouin zone" (of height Hlo) 
with states originating from the lowest three Bloch bands 
for w/(27r) — 5.30 kHz, the frequency of the extraordi- 
nary peak in Fig. and — 0.7, 1.0, and 1.3. There 
are various avoided crossings indicating multiphoton-like 
couplings between the bands; however, with Afc = 0.1 fcL 
the wave packet evolving from the initial distribution (jS]) 
mainly explores the interval of quasimomenta indicated 
by the shaded areas. The quasienergy band originating 
from the lowest unperturbed energy band n = 1 is shown 
enlarged in the insets; with increasing Kq this band is 
pierced through from below by the quasienergy band 
n = 2, displaced down in energy by ifuo against that rep- 
resentative which is continuously connected to the bare 
n = 2 Bloch band. This penetration results in "active" 
avoided crossings signaling a strong-field-induced three- 
photon resonance; this is responsible for the anomalous 
peak at w/(27r) = 5.30 kHz. 

The dynamics underlying that peak should thus be dis- 
cussed in terms of the morphology of the surfaces which 
emerge when the quasienergies are considered as func- 
tions of both the wave number k and the instantaneous 
amplitude F (or Kq): When the driving amplitude F{t) 
increases during the upward ramp of a pulse, the ini- 
tial distribution is shifted almost adiabatically on its 
quasienergy surface, parallel to the Kq axis. As long as 
the maximum value of lies below the critical regime 
where this surface is first being pierced by another one, 
the initial distribution is restored with only minor dis- 
tortion when the amplitude returns to zero, resulting in 
negligible escape probability. However, when the moving 
distribution hits an avoided-crossing regime, part of the 
wave function undergoes a Landau-Zener-type transition 
to the anticrossing band. Both parts of the wave func- 
tion then evolve separately on their respective surfaces, 
until they meet for a second time during the downward 
ramp, when they interfere and thereby establish the fi- 
nal occupation probabilities of the bands involved. This 
mechanism of splitting and interference implies that there 
should be Stiickelberg-like oscillations when the final oc- 
cupation probabilities are monitored while the length of 
the pulses' plateau segment is varied, because varying 
the plateau duration means varying the relative phase 
picked up by the two interfering components. Indeed, 
these oscillations are clearly visible in Fig. [S) 

We remark that the standard perturbative m-photon 
resonances can be grasped in a similar manner: For fre- 
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FIG. 5: (Color online) Stiickelberg oscillations of the escape 
probability in response to prolongation of the plateau dura- 
tion thoid, for oj/{2-k) = 5.30 kHz, and isT""''' = 1.0 (dashed) 
and 1.3 (solid line). 

quencies such that Eq. ([9]) holds, two quasienergy sur- 
faces are degenerate already at F = 0, so that adia- 
baticity is disabled and the wave function splits right 
at the beginning of a pulse [l^ . Seen against this back- 
ground, a perturbative resonance corresponds to the re- 
moval of a quasienergy degeneracy already present at 
= 0, while a nonperturbative one emerges when ac- 
Stark-shifted Bloch bands penetrate each other at a cer- 
tain finite driving strength. 



IV. CONCLUSIONS 

When viewing a time-periodically forced optical lattice 
as a spatiotemporal crystal, the natural basis states are 
the spatiotemporal Bloch waves ([TOl) : the energy bands 
En{k) of the undriven lattice turn into quasienergy bands 
e„(fc). The latter depend not only on the lattice param- 
eters, but also on the parameters of the driving force. 
While they differ barely from the unperturbed energy 
bands as long as the driving amplitude is weak, corre- 
sponding to values Kq ^ 1 of the dimensionless para- 
meter (j4|), they become strongly distorted, and even pen- 
etrate each other, in the nonperturbative regime. 

When subjected to pulsed forcing with an amplitude 
which changes slowly compared to the period T = 2-k /lo 
of the drive, a wave packet can adjust itself adiabatically 
to a mere distortion of its quasienergy band. However, 
when the wave packet explores a part of a quasienergy 
band which is pierced by another one, as exemplified in 
Fig. |4l Landau-Zener transitions occur; this mechanism 
leads to strong nonperturbative resonances at frequencies 
not given by the simple condition In principle, such 
resonances should also occur in solids irradiated by strong 
laser pulses; however, there they would be masked by a 
host of competing effects. The experimentally proven 
good controllability of ultracold atoms in forced optical 
lattices makes such systems a far better testing ground 
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for these dynamics. 

Our study has been restricted to the single-particle 
level; it is reasonable to expect that the phenomena ex- 
emplarily discussed in the present work can immediately 
be detected with sufficiently dilute or close-to-ideal Bose- 
Einstein condensates in driven optical lattices [12]. Even 
more, it appears equally feasible to perform the exper- 
iments suggested here under conditions of sizable inter- 
particle interactions, or even of strong correlations. The 
question how the single-particle scenario outlined above 
is modified then opens up far-reaching further lines of 



investigation, concerning both experiment and theory. 
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